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Abstract 

In the first paper of this series [S. Torquato, J. Chem. Phys. 136, 054106 (2012)], analyti- 
cal results concerning the continuum percolation of overlapping hypcrparticles in d-dimensional 
Euclidean space M*^ were obtained, including lower bounds on the percolation threshold. In the 
present investigation, we provide additional analytical results for certain cluster statistics, such as 
the concentration of A;-mers and related quantities, and obtain an upper bound on the percolation 
threshold r]c- We utilize the tightest lower bound obtained in the first paper to formulate an effi- 
cient simulation method, called the rescaled-particle algorithm, to estimate continuum percolation 
properties across many space dimensions with heretofore unattained accuracy. This simulation 
procedure is applied to compute the threshold r]c and associated mean number of overlaps per 
particle Mc for both overlapping hyperspheres and oriented hypercubes for 3 < d < 11. These 
simulations results are compared to corresponding upper and lower bounds on these percolation 
properties. We find that the bounds converge to one another as the space dimension increases, but 
the lower bound provides an excellent estimate of rjc and Mc, even for relatively low dimensions. 
We confirm a prediction of the first paper in this series that low-dimensional percolation properties 
encode high-dimensional information. We also show that the concentration of monomers dominate 
over concentration values for higher-order clusters (dimers, trimers, etc.) as the space dimension 
becomes large. Finally, we provide accurate analytical estimates of the pair connectedness function 
and blocking function at their contact values for any d as a function of density. 
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I. INTRODUCTION 



In the first paper (paper I)i of this series of two papers, we obtained a number of ana- 
lytical results concerning the continuum percolation of overlapping hyperspheres and over- 
lapping oriented hypercubes that applied across all Euclidean space dimensions. Among 
other results, it was shown analytically that certain lower-order Fade approximants on the 
mean cluster number S are lower bounds on S for both hyperspheres and hypercubes in 
d- dimensional Euclidean space M*^ and that they become exact asymptotically as (i — )■ oo. 
In this asymptotic limit, the dimensionless density at percolation rj^ tends to 2~'^. An im- 
portant consequence of the aforementioned analysis is that this large-c? percolation value 
is an important contribution to the low- dimensional percolation value. In other words, 
low- dimensional results encode high- dimensional information. Percus-Yevick-like approxi- 
mations for the cluster number S were also obtained that also become asymptotically exact 
as d — 7- oo. The analysis was aided by a striking duality between the equilibrium hard- 
hypersphere (hypercube) fluid system and the continuum percolation models of overlapping 
hyperspheres (hypercubes), namely, 

P{r;ri) = —h{r; —1]) (1) 

where P{r; t]) is the pair connectedness function at some radial distance r and reduced 
density rj for the continuum percolation models and h{r\ri) is the total correlation function 
for the equilibrium hard-particle models.- It was shown that the large-c? percolation threshold 
rjc of overlapping hyperspheres and hypercubes is directly related to the large-c? freezing-point 
density of corresponding equilibrium hard-particle models.^ The extension of these results 
for overlapping hyperspheres and hypercubes to the case of overlapping particles of general 
anisotropic shape in d dimensions with a specified orientational probability distribution was 
also described. 

The bounds and approximations reported in paper I were applied to assess the accuracy 
of previous computer simulation results for r/c that span dimensions up to = 20 in the case 
of overlapping hyperspheres^!^ and up to c? = 15 in the case of hypercubes.- It is convenient 
to restate the best lower bound on r/^ that was found in Ref. namely. 



where C3 and C4 are the trimer and tetramer statistics defined and computed as a function 
of d in paper I. Comparison of this lower bound to Kriiger's simulation data for hyperspheres 
revealed that the bound became progressively tighter as d increased and became an excel- 
lent estimate for d > 7. Since it becomes increasingly challenging to estimate percolation 
thresholds from simulations in high dimensions, it was not surprising that even Kriiger's 
high-quality simulation data fell slightly below the lower bound (|2]) for 8 < c? < 11. The 
analytical results of paper I revealed that the simulation data reported in Ref. |4| for both hy- 
perspheres and hypercubes were considerably more problematic. In particular, the authors 
reported incorrectly that the quantity for these systems were nonmonotonic in dimen- 
sion and that hyperspheres have lower thresholds than hypercubes in higher dimensions 
while the reverse is true in lower dimensions. 

One of the main purposes of the present paper is to exploit the accuracy of the lower 
bound ([2]) to provide an efficient simulation method, called the rescaled-particle algorithm, 
to estimate continuum percolation thresholds across many dimensions. Another objective is 
to obtain an upper bound on rjc as well as to provide additional analytical results for certain 
cluster statistics. In Sec. [Tll we derive these analytical results. In Sec. IIIIl we describe in 
detail the rescaled-particle method, which is applicable for general continuum percolation 
models (e.g., spherical and nonspherical particle shapes). Results for the thresholds of both 
overlapping hyperspheres and oriented hypercubes for 3 < (i < 11 are compared to upper 
and lower bounds on rjc in Sec. IIVI In Sec. |Vl we summarize our conclusions and discuss 
future work. 

II. ADDITIONAL ANALYTICAL RESULTS 
A. Basic Definitions 

A prototypical continuum percolation model consists of equal-sized overlapping (Poisson 
distributed) hyperparticles in R'^ at number density p; see paper I and numerous references 
therein. It is convenient to introduce the reduced number density 77, defined by the relation 

V = pvi, (3) 

where Vi is the 0?- dimensional volume of a hyperparticle; see Ref. Q for an explicit expression 
of this quantity for a hypersphere, for example. A cluster statistic that has been considered 
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by various investigators in one, two and three dimensions for overlapping spheres is n^, the 
average number of /c-mers per unit number of particles.-""- This k-mer statistic obeys the 
following constraint 

^knk = l, (4) 

k=l 

where it is to be noted that pk = kuk is the probability that a given particle is part of a 
k-meiA In what follows, we will derive estimates for Uk and related cluster statistics for 
arbitrary dimension. 

It was shown in Ref. ^ that can be explicitly expressed for any d for overlapping 
hyperspheres as certain multidimensional integrals involving exponentials whose arguments 
contain the union volume of k "exclusion spheres." Two spheres of radius D/2 are considered 
to be connected if they overlap, i.e., if the center of one lies within a spherical "exclusion" 
region of radius D centered around the other sphere (see Fig. 3 of paper I). For example, 
for k = 1, k = 2 and A; = 3, we have 

ni = exp[-2'^ri] (5) 

P 
2 



n2 = ^ I exp[-pv2{r; D)]f{r)dr (6) 



^3 = ^ / ^^"12 / c?ri3 /(ri2)/(r23)exp[-pt;3(ri2,ri3,r23;L')] 

+ 1" / '^^12 / rfri3exp[-pt;3(ri2,ri3,r23;L')], (7) 

•J Jr^ JRd 

where f„(ri, r2, . . . , r„; D) is the union volume of n spheres of radius D centered at positions 
ri, r2, . . . , r„, Tij = rj — rj and Vij = \rij\ {i ^ j). Moreover, the radial function /(r) defines 
the connectedness criterion, i.e., 

f{r) = Q{D-r), (8) 

and 

is the Heaviside step function. Note that the factor 2°^, appearing in Eq. ([5]), is the ratio of 
the exclusion volume Vex to the volume of a sphere. By virtue of the fact that the spheres 
are Poisson distributed in space, it follows that the mean number of overlaps per sphere J\f 
is given by 

Af = pvex = (10) 
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The average number of clusters per unit volume pc is directly related to the sum over n^,- 
namely 

oo 

7 = E-- (11) 

^ k=l 

Note that the ratio pc/ p is the average number of clusters per particle and its derivative 
with respect to p (or rj) determines the contact value of the blocking function B{r),— which is 
directly related to the conditional probability of finding two particles belonging to different 
clusters separated by distance r, given that one of the particles is at the origin. For example, 
in the case of hyperspheres of diameter D, the three-dimensional expression given in Ref. l9 
generalizes as follows: 

rl( n I n\ , , 

= ~2'^-^B{D). (12) 



drj 

We recall here that for overlapping hyperspheres, the blocking function B{r) can be obtained 
immediately from the pair connectedness function P(r) for any radial distance r via the 
relation^ 

P{r) + B{r) = 1 . (13) 

The fact that P(r) is bounded in the interval [0, 1], implies the same bounds on B{r). 

The average cluster number Q is the average number of particles in a randomly chosen 
cluster and is the inverse of pc/ Pr namely, 

-1 



oo 



^ = -=11.^^1 ■ (14) 

This is to be distinguished from the cluster number 5* (average number of particles in the 
cluster containing a randomly chosen particle), which is related to the second moment of 



S = '^k^nk, r]<r]c. (15) 



k=l 



Unlike 5, which can also be expressed in terms of the pair connectedness function P{r),— 
the average cluster number Q does not diverge at the percolation threshold when d > 2.- 
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B. Estimates of Cluster Statistics and Upper Bounds on rjc and Mc 

It was shown in Ref. Q that the mean number of overlaps per sphere at the threshold Afc 
[cf. ffTU]) ] tends to unity as d — )■ oo, i.e., 

A/'c = ''7c— ~ 1, d — )■ oo, (16) 

Vl 

which applies to spherical as well as nonspherical particles (with specified orientational 
distribution). This does not mean that the concentrations of monomers, dimers, trimers, 
etc. at the threshold are negligibly small, even if finite clusters become more ramified as the 
space dimension grows. ^ To explicitly prove this property, we first observe that the exact 
formula for rii [cf. ([5])] together with the exact asymptotic result 

Vc^^, d^oo, (17) 

which applies to any oriented centrally symmetric particle (e.g., spheres, cubes, ellipsoids, 
etc.)i, imphes the following asymptotic result 

ni ~ exp(-l) = 0.3678794. . ., d ^ oo. (18) 

We will now show that as d becomes large, monomers, as opposed to any fc-mer for k > 2, 
are dominant in so far as concentration is concerned, i.e., rii is appreciably larger than n2 
and therefore is appreciably larger than with A; > 3, since Uk > Uk+i for any positive 
but bounded rj. Let us begin with the formula ([6]) for the dimer statistic n2, which we can 
rewrite as follows: 

where 

- ~ Twwrm Jo "^^^ ^ 

Here we have used the fact that V2{r; R) = 2vi{R) — v^^{r;R), where the latter quantity 
is the intersection volume of two spheres of radius R whose centers are separated by the 
distance r. The dimensionless intersection volume a;(r; R), which has support in the interval 
[0, 2R], has been explicitly given for any d in a variety of representationsi*^ and played an 
important role in paper I. Now since a{D; D) decays to zero exponentially fast according to 



d2'^~^r]exp[-2'^+^7]] [ r"'-^ exp[2 V(^; ^)], (19) 
Jo 

-int(r;/2) 21(1 + ^/2) r-os-Hr/i2R) 
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the asymptotic relation^ 



„(C;C)~(«V'7^)'"-lj. (21) 



TT 



it immediately follows from ( IT9l) and (ITTl) that 



and hence we find 



^ ^^P^ = 0.06766764 . . . , d ^ oo, (22) 



— ~2exp(l) = 5.436563..., ^ oo, (23) 



^2 

which is what we set out to prove. 

We now derive lower bounds on for A; > 2 as a function of d for any rj. Let us begin 
with the case k = 2. Since a{r; D) is a monotonically decreasing function of r, we have that 
a{D; D) < a{r; D) in the interval [0,D] and hence combined with the exact formula (fT9|l 
yields the lower bound 

n2 > 2^-^77 exp[-2"'+^r/] exp[2V(^; D)]. (24) 



It is noteworthy that in light of ( ]2T1) . the lower bound ( ]24l) becomes asymptotically exact 
in the high-d limit, i.e., we recover (!22|) . Using similar arguments and the formulas for 
given in Ref. ^, we obtain the following generally weaker lower bounds on rij. for any k: 

2{k-l)d 

rik > ^j^_^y V'~' exp[-fc2^r/]. (25) 



While for k = 1, this bound is exact, it is weaker than fl24l) for k = 2. 

Through second order in t], formulas ([5]), and ([7]) for rii, n2 and n^, respectively, yield 

m = 1 - 2'^7] + 2^'^'^r]'^ + 0{r]^) (26) 
n2 = 2''-'r] - (^2^^ + + ^(^') (27) 

)2d-l I ^3 A 2 I /'O/'^Sn 



n3 = (^22'^-i + yJ ^ + ) (28) 

Similarly, using expression (fT4|) and the relations immediately above, we can obtain the 
corresponding density expansion for Q: 



(29) 
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It was noted in Ref. Qthat the pole of the [1,1] Fade approximant of the density expansion 
of Q for (i = 3 yielded an upper bound on the threshold rjc for overlapping spheres. For 
general d, the [1,1] Fade approximant for Q for either overlapping hyperspheres or oriented 
hypercubes is given by 

where the trimer statistic C3 for both models is given in Ref. Ill We now observe that the 
pole of (|5U]) is an upper bound on rjc for overlapping hyperspheres and oriented hypercubes 
for any (i > 3, i.e., 

- ri 1^ c, y d>3. (31) 

L-L -r e.22(d-i) J 

Using the same methods described in paper I, it is straightforward to prove that for suf- 
ficiently large d and any rj < rjc, relation ( 130|) bounds Q from below and hence (13T!) is a 
rigorous upper bound on the threshold. This is consistent with the observation (noted in 
Sec. Ill All that the actual function Q does not diverge when the mean cluster number 5" 
diverges, i.e., when r] — )■ rjc- We will show in Sec. [IV] that the expression fl^ bounds 
the simulation data for the threshold from above for both hyperspheres and hypercubes 
for 3 < d < 11. Moreover, according to Ref. Q, the second term within the brackets of 
inequality f l3T|) goes to zero exponentially fast in the limit d — )■ 00, and hence this bound 
asymptotically becomes 

Vc<^, d^oo, (32) 

which is the exact asymptotic result.- Since the lower bound ([2]) also becomes exact in this 
high-ci limit, the bounds ([2]) and ( l3T]) converge to the exact asymptotic value of 2"*^. 

Finally, we note that combining relation (fT4|) and the approximant ( 130|) gives the following 
approximation of Pc/ p: 



1 - [2'^"! + V 



r . (33) 

Substituting ( !33l) into ( fT2l) yields the following the approximation for the contact value of 
the blocking function 

B{D) = 1 - P{D) ^ ^ . (34) 

1-^r/ 

In the Appendix, we provide plots of B{D) and F(D) versus r] for selected dimensions. 
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III. EFFICIENT ALGORITHM TO COMPUTE r/c ACROSS ALL DIMENSIONS 



A. Particle- Addition Method 

A commonly employed approach to estimate the percolation threshold rjc for continuum 
percolation in two and three dimensions is the particle-addition method.—"— Starting from a 
configuration of a small number of particles with random positions in the simulation domain 
subject to periodic boundary conditions, new particles are added to the domain sequentially 
with randomly chosen positions in the simulation domain. Each time a new particle is added, 
the largest cluster in the system (i.e., the one containing the largest number of particles) 
is identified using a burning-algorithm.— This process is repeated until a system-spanning 
cluster forms. 

Although conceptually intuitive and used across dimensions,-"^ this method become pro- 
gressively less computationally efficient as the dimension increases.- First, to obtain accurate 
estimates of the percolation threshold, the size of the particles should be much smaller than 
the linear extent of simulation domain so that adding a single particle leads to a very small 
increase in the reduced density rj. Without a priori knowledge of the percolation thresh- 
old, one needs to start with sufficiently dilute particle configurations, i.e., with very small 
77. Therefore, an extremely large number of particles needs to be added to the simulation 
domain until the system percolates. For each particle addition, one needs to identify the 
largest cluster and check whether it spans the system, which makes the method computa- 
tionally very expensive. Moreover, the dynamic nature of particle addition makes it difficult 
to implement efficient methods to check for local particle connectivity (overlaps of pairs 
of particles, e.g., the cell method). This problem increases in severity as the dimension 
increases. 

We note that highly efficient algorithms have been developed for investigating clustering 
and percolation properties of overlapping disks in two dimensions. For example, several 
frontier-tracking methods have been devised to provide very precise estimates of the per- 
colation threshold rjc of overlapping disks.—"— However, such algorithms cannot be applied 
in higher dimensions because there is no analogous localized boundary in the gradient per- 
colation method for high- dimensional systems.— Other efficient variations of the particle 
addition method have been implemented, e.g., particles are added only to a single growing 
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cluster, and the cluster-size distribution rather than the spanning cluster is used to obtain 
rjc'r^ and rapid "union-find" methods have been developed to keep track of the connected 
clusters as particles are added.— However, these variations of the particle-addition method, 
although much more efficient than the original one, becomes progressively more difficult to 
apply as the dimension increases. 

B. Rescaled-Particle Method 

Tight rigorous lower bounds on the percolation threshold, such as the ones derived in Ref. 
1, enable us to devise a highly efficient method to estimate rjc for overlapping particles with 
arbitrary shapes and orientations in M.'^. The basic idea is to generate initial static particle 
configurations at a value of r] that is taken to be the best lower-bound value, allowing one 
to finally arrive at the critical value in a computationally efficient manner, even in high 
dimensions. We note that procedures in which particles are rescaled to obtain rjc have been 
previously proposed.— However, to the best of our knowledge, our method is the first one to 
combine the particle rescaling procedure with the tightest lower bound values to efficiently 
and accurately obtain rjc- 

Initially, a Poisson distribution of a large number of points in the simulation domain 
is generated. Each point in configuration is then taken to be the centroid of a particle 
with a specified or random orientation and a characteristic particle length scale io (e.g., 
the diameter of a sphere). The initial value of io is chosen such that the reduced density 
1] of the system equals the tightest lower bound value. For each particle i, a near-neighbor 
list (NNL) is obtained that contains the centroids of the particles j whose distance Dij to 
particle i is smaller than 7^0 (7 > !)• The value of 7 generally depends on the continuum- 
percolation models of interest. A rule of thumb for a choosing good value of 7 is that the 
NNL list only contains particles that overlap at the percolation threshold. In our simulation 
for hyperspheres and hypercubes, we have used 7 G [1.05,1.5], depending on the space 
dimension. Then the particle sizes are slowly and uniformly increased by increasing io (i.e., 
rescaling the particles), which leads to an increase of the reduced density by a small amount 
6t]. After each rescaling, the particles in the NNL are checked for overlap and the largest 
cluster in the system is identified. The process is repeated until a system-spanning cluster 
forms. 
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FIG. 1: Cluster statistics, including the linear size of the largest cluster— (a) and the fraction of 
particles in the largest cluster (b), associated with the particle-addition and rescaled-particle meth- 
ods for overlapping hyperspheres in M^. The different methods clearly produce very similar cluster 
statistics yet the rescaled-particle method is much more computationally efficient, as explained in 
the text. 

Since static particle configurations and predetermined NNL are used, the complexity for 
identifying clusters is significantly reduced. In addition, as the space dimension increases, 
it was shown that the lower bounds derived in paper I become increasingly tighter.- Thus, 
total number of rescaling before percolation is achieved is much smaller than the total 
number of particle additions, which dramatically improves the efficiency of the algorithm. 
Furthermore, the increase of the reduced density 61] can take arbitrarily small values when 
7]c is approached, rather than a fixed discontinuous value determined by the particle size 
via particle-addition methods. This smooth approach to the critical value allows a more 
accurate estimate of the percolation threshold. 

Similarly, given an upper bound on rjc, such as inequality ( 13T]) . one can start with a 
percolated system, and rescale the particles (i.e., decreasing the particle size Eq) to reduce 
7]. Figure [1] shows the cluster statistics (e.g., the linear size of the largest cluster^^ and the 
fraction of particles in the largest cluster) associated with the particle-addition method and 
the rescaled-particle method for overlapping hyperspheres in M'^. The initial configurations 
for the rescaled-particle method include both a non-percolated configuration with t] equal to 
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the lower-bound value riL and a percolated configuration with t] equal to the upper-bound 
value rjij. We note that close to percolation, the cluster containing the largest number of 
particles also possesses the largest linear size. 

It is clear that the particle-addition method and the rescaled-particle method produce 
very similar cluster statistics yet the rescaled-particle method is considerably more computa- 
tionally efficient in three dimensions. The computational efficiency improves for overlapping 
hyperspheres and hypercubes as the space dimension increases beyond three. Since the up- 
per bound on rjc derived here is not as tight as the lower bound in relatively low dimensions, 
we mainly use the rescaled-particle method starting from non-percolated configurations with 
7] equal to the tightest lower-bound value, as predicted from ([2]). However, as noted in Sec. 
IIIBl as d — )■ oo, both the upper and lower bounds converge to the exact asymptotic value 
of 2-'^. 

IV. SIMULATION RESULTS FOR OVERLAPPING HYPERSPHERES AND 
ORIENTED HYPERCUBES 

Using the rescaled-particle method starting from a reduced density given by the lower- 
bound estimate ([2]), we compute the percolation threshold for overlapping hyperspheres 
and oriented hypercubes in dimensions two through eleven. For each dimension, different 
system sizes N are used. Specifically, we employ N = 10000, 50000, 100000 for d = 2,3,4, 
N = 50000, 100000, 500000 for d = 5,6,7, N = 100000, 500000, 1000000 for d = 8,9, 10, 
and = 1000000, 2500000, 5000000 for c/ = 11 and the results are extrapolated by spline 
fitting the finite-system-size data in a log-log plot to obtain the infinite-system-size estimate 
of rjc- For each system size, the percolation threshold is obtained by averaging over 1000 
independent particle configurations for d = 2 and 3, 500 independent particle configurations 
for 4: < d < 8, and 100 independent particle configurations for 9 < < 11. 

The obtained percolation threshold values rjc for overlapping hyperspheres and hypercubes 
in dimensions two through eleven are respectively given in Table 1 and Table 2, and displayed 
in Figs. |2] and [3l We also provide in the tables and figures the corresponding values of the 
lower bound (E]) and upper bound ( I3T1) on t]c for purposes of comparison. Note that our 
simulation data lie very close to the lower-bound values, and that the lower bounds and 
data converge quickly to one another as d increases. Moreover, we include the numerical 
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TABLE I: Estimates of the percolation threshold ij^. for overlapping hyperspheres as obtained from 
the rescaled-particle algorithm, the lower bound ([2|), and the upper bound ([3T]) . Also included 
are the numerical estimates rj* of the percolation threshold from a previous study'^ that satisfy 
the bounds ([2]) and (I3ip . These results are not reported for d > 8 because they violate the lower 
bounds. 



I 

a 


Vc 


Vc 


Vc 


Tc 


2 


0.748742. . . 


1.1282 


1.12810(3) 




3 


0.271206. . . 


0.3418 


0.34289(2) 


0.363636. . . 


4 


0.111527. .. 


0.1300 


0.1304(5) 


0.167373. . . 


5 


0.0488542. . . 


0.0543 


0.05443(7) 


0.0788179. . . 


6 


0.0222117. .. 


0.02346 


0.02339(5) 


0.0376720. . . 


7 


0.0103452. . . 


0.0105 


0.01051(3) 


0.0181921. .. 


8 


0.00489917. . . 




0.004904(6) 


0.00885075. . . 


9 


0.00234800. . . 




0.002353(4) 


0.00432995. . . 


10 


0.00113534. . . 




0.001138(3) 


0.00212726. . . 


11 


0.000552682. . . 




0.0005530(3) 


0.00104854. . . 



estimates of ri^ from previous simulation studies for hyperspheres^ and hypercubes'*^ in case 
they do not violate the bounds ([2]) and ( ISTl) . It can be clearly seen that our rescaled-particle 
method yields much more accurate estimates of t^c, especially in high dimensions, since we 
start with particle configurations that are already very close to percolation. 

Following Ref. Q we use the threshold estimate obtained from the [2, 1] Fade approximant 
of S as the basis to obtain accurate analytical approximations for r^c that applies across all 
dimensions for hyperspheres and oriented hypercubes. Specifically, we fit the following 
function to the simulation data for 2 < < 11: 



where 



(2) 
^0 



6^ 62 



C3 
1 + — 



(2) 



2Cq C A 
1 H H 



(35) 



(36) 
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TABLE II: Estimates of the percolation threshold r/^ for overlapping hypercubes as obtained from 
the rescaled-particle algorithm, the lower bound ([2|), and the upper bound (j3ip . Also included 
are the numerical estimates rj* of the percolation threshold from a previous study^ that satisfy 
the bounds ([2]) and (I3ip . These results are not reported for d > 5 because they violate the lower 
bounds. 



2 


0.732558. . . 


1.098 


1.0982(3) 




3 


0.256680. . . 


0.3248 


0.3247(3) 


0.347824. . . 


4 


0.103286. .. 


0.12 


0.1201(6) 


0.158416. .. 


5 


0.0447161. . . 




0.05024(7) 


0.0742456. . . 


6 


0.0202386. . . 




0.02104(8) 


0.0354571. . . 


7 


0.0094301. . . 




0.01004(5) 


0.0171512. .. 


8 


0.00448213. . . 




0.004498(5) 


0.00837119. . . 


9 


0.00216025. . . 




0.002166(4) 


0.00411207. . . 


10 


0.00105159. . . 




0.001058(4) 


0.00202930. . . 


11 


0.000515602. . . 




0.0005160(3) 


0.00100485. . . 



0.2 
0.15 

< 

0.1 
0.05 



"4 5 6 7 8 9 10 n 
d 

FIG. 2: Percolation threshold rjc versus dimension d for overlapping hyperspheres as obtained from 
the lower bound ([2]), the upper bound (f3T]) and the simulation data. 
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FIG. 3: Percolation threshold 7]c versus dimension d for overlapping hypercubes as obtained from 
the lower bound ([2]), the upper bound ([3T]) and the simulation data. 

is the pole associated with the [2, 1] Pade approximant for the mean cluster numver S, 
i.e., the tightest lower bound for rjc, explicitly given by Eq. (119) of paper I. This lower 
bound becomes exact for sufficiently large d. Therefore, in agreement with the conclusions of 
Ref. [l|, we see that the high-dimensional percolation behavior is an important contribution 
to the low-dimensional percolation value. In other words, low- dimensional results encode 
high- dimensional information. 

We find that bi = 2.45074 and 62 = —1.65036 for hyperspheres with correlation coeffi- 
cient equal to 0.993194; and bi = 2.57917 and 62 = —2.29755 for hypercubes with correlation 
coefficient equal to 0.992262. Note that due to the quality of the available numerical data 
reported in paper I, the analogous analytical approximations based on those data for hy- 
perspheres and hypercubes only used numerical threshold estimates for 2 < d < 7 and 
2 < d < 4, respectively. Thus, the formula (l35i) for hyperspheres and hypercubes supersedes 
in accuracy the ones provided in Ref. Q. 

As pointed out in Ref. Q, the numerical estimates of rjc for hypercubes in 5 < < 15 by 
Wagner et alA violate the tightest lower bound ([2]). In addition, these simulation data are 
questionable in high dimensions since the mean number of overlaps per particle M [defined 
in Eq. (fTOl) ] evaluated at the percolation threshold, i.e., Mc = 2'^rjc is incorrectly found to 
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FIG. 4: The mean number of overlaps per hypersphere at percolation threshold Nc = '^'^f]c as a 
function of d. Also shown are the quantities 2'^r]L and 2'^riu, where rji and rju are respectively the 
tightest lower bound ([2]) and upper bound ([3T]) . It is clear that Mc is a monotonic function of d 
and quickly converges to the asymptotic value of unity as d increases. 

be a nonmonotonic function of d. In particular, these authors found that A/'c first decreases 
as d increases for 2 < < 9 and then increases as d increases for 10 < d < 15. This led to 
the incorrect conclusion that hyperspheres have lower thresholds than hypercubes in higher 
dimensions while the reverse is true in lower dimensions. 

In Figs, m and El we show Afc as a function of d computed using the estimates of r]c 
obtained from our simulations for overlapping hyperspheres and hypercubes, respectively. 
The quantities and are also shown for purposes of comparison, where rj^ and rju 
are respectively the tightest lower bound ([2]) and upper bound (!3T|) . It can be clearly seen 
that A/'c for both overlapping hyperspheres and hypercubes are indeed monotonic functions 
of d, which quickly converge to the asymptotic value of unity as d increases. This indicates 
again that the large-d asymptotic percolation value is an important contribution to the low- 
dimensional percolation value. Moreover, one can see from Tables 1 and 2 that hypercubes 
always have a lower threshold than hyperspheres for any fixed finite dimension, and the 
thresholds of these two systems approach one another in the limit — )■ oo, as predicted in 
Ref. ll. 

In Fig. [6], we plot the concentration of monomers and dimers, ni and n2, at the percolation 
threshold 7]c as a function of dimension d for overlapping hyperspheres as obtained from the 
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FIG. 5: The mean number of overlaps per hypercube at percolation threshold Afc = as a 
function of d. Also shown are the quantities 2'^r]L and 2'^riu, where rji and rju are respectively the 
tightest lower bound and upper bound (f3T]) . It is clear that Mc is a monotonic function of d 
and quickly converges to the asymptotic value of unity as d increases. 

exact expressions (|5]) and (fT9|) and the simulation data for rjc given in Table I. Observe that, 
consistent with analysis given in Sec. Ill Bt rii becomes appreciably larger than n2 (and hence 
ns,n4, etc.) as the dimension increases. We include in the figure the lower bound (|24l) on 
n2, which we see becomes tighter as the space dimension increases, as the analysis of Sec. 
Ill Bl predicts. 

V. CONCLUSIONS AND FUTURE WORK 

We have supplemented the analytical results obtained in paper I by deriving additional 
formulas and bounds for certain cluster statistics, such as the concentration of fc-mers and 
related quantities, and obtained an upper bound on the percolation threshold rjc- We utilized 
this upper bound and the tightest lower bound on rjc obtained in paper I to devise an efficient 
simulation method, called the rescaled-partide algorithm, to estimate continuum percolation 
properties across many space dimensions. We applied this simulation procedure here to 
compute, with heretofore unattained accuracy, the threshold rjc and associated mean number 
of overlaps per particle Afc for both overlapping hyperspheres and oriented hypercubes for 
3 < d < 11. Comparison of these simulations results to corresponding upper and lower 
bounds on these percolation properties revealed that the bounds converge to one another as 
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FIG. 6: Monomer and dimer concentrations, ni and n2, at the percolation threshold r/c as a function 
of dimension d for overlapping hyperspheres as obtained from the exact expressions ([5]) and (jl9p 
and the simulation data for tJc given in Table I. Included in the figure is the lower bound (j'24|) on 
n2. 

the space dimension increases. It is noteworthy that the lower bound provides an excellent 
estimate of rj^ and A/'c, even for relatively low dimensions. We confirmed a prediction of 
paper I that low- dimensional percolation properties encode high- dimensional information. 
We also showed that the concentration of monomers dominate over concentration values for 
higher-order clusters (dimers, trimers, etc.) as the space dimension becomes large. Finally, 
we provided accurate analytical estimates of the pair connectedness function and blocking 
function at their contact values for any d as a function of density. 

In paper I, the extension of the continuum percolation results obtained for overlapping 
hyperspheres and oriented hypercubes to cases in which the overlapping hyperparticles are 
nonspherical (anisotropic in shape) with some specified orientation distribution function 
(e.g., random orientations) was briefly discussed. Future work will expound on this exten- 
sion to overlapping anisotropically-shaped hyperparticles with random orientations.. The 
exploration of the generalizations of the techniques of this series of papers to bound per- 
colation thresholds in the lattice setting^' (bond and site percolation^'^>2i) represents an 
intriguing area for future research. 
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Appendix: Approximations for the Blocking and Pair Connected- 
ness Functions at Contact 

In Sec. II, we noted that the inverse of the cluster number Q is the of the average number 
of clusters per unit volume pc/p [see Eq. f[T^ ]. We derived the approximation ( 15U]) for Q, 
which immediately leads to the approximation ( 133|) for Pc/p and thus, the approximation 
for the contact values of the blocking function B{D) and pair connectedness function 
P{D). 

In this appendix, we compare the approximation ([3SD for pc/ p versus rj with available 
numerical data for overlapping spheres in M^.- These results are plotted in Fig. [71 Observe 
that the approximation agrees with the simulation data very well and bounds the data from 
below. 
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FIG. 7: The approximation (I33p of Pc/p as a function of r? for three-dimensional overlapping 
spheres. Also shown are the simulation data reported in Ref. |9| for purposes of comparison. It can 
be seen that the approximation agrees very well with the simulation data and bounds the data 
from below. 



Since the approximation fl33|) for pc/ p improves as the space dimension increases beyond 
three, we expect that results derived from it, such as the relation for the blocking 
function and pair connectedness function at contact, will provide accurate approximations 
across dimensions for d > A. In particular, we provide plots of the contact values B{D) 
and P{D) = 1 — B{D) [as predicted by flM|) ] versus r] for overlapping hyperspheres in 
dimensions 3, 7 and 11 up to the respective percolation thresholds in Figs. [H [9] and [TOl The 
functions B{D) and P{D) are equal to unity and zero, respectively, at r] = and decrease 
and increase monotonically with increasing rj up to 1]^. We also see that B{D) and P{D) 
vary less appreciably with increasing rj from the maximum value of unity and minimum 
value of zero, respectively, as the space dimension increases. 
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FIG. 8: The contact values of the blocking function B{D) and the pair-connectedness function P{D) 
versus i] up to ?7c = 0.34289 (see Table I) for overlapping spheres in . Note that B{D)+P{D) = 1 

[c.f. m]- 
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FIG. 9: The contact values of the blocking function B{D) and the pair-connectedness function 
P{D) versus rj up to rjc = 0.01051 (see Table I) for overlapping hyperspheres in M^. Note that 
B{D) + P{D) = 1 [c.f. dMl)]. 
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FIG. 10: The contact values of the blocking function B{D) and the pair-connectedness function 
P{D) versus r/ up to r/c = 0.0005530 (see Table I) for overlapping hyperspheres in M^^. Note that 
B{D) + P{D) = 1 [c.f. (IMI)]. 
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